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Wide-angle, non-imaging air Cherenkov detectors provide a way to observe cosmic gamma-rays which is complementary to obser- 
vations by imaging Cherenkov telescopes. Their particular strength lies in the multi-TeV to ultra high energy range (Ey > 30 TeV), 
where large effective areas, yet small light sensitive areas per detector station are needed. To exploit this potential to full extent, a 
large station spacing is required to achieve a large effective area at a reasonable effort. 

' In such a detector, the low number of signals per event, the absence of imaging information, and the poor signal to noise ratio 
of Cherenkov light to night sky brightness pose considerable challenges for the event reconstruction, especially the gamma hadron 
separation. The event reconstruction presented in this paper has been developed for the wide-angle detector HiSCORE, but the 
concepts may be applied more generically. It is tested on simulated data in the 10 TeV to 5 PeV energy range using the air shower 
simulation CORSIKA and the HiSCORE detector simulation sim_score. For the tests, a regular grid of 22 x 22 detector stations 
with a spacing of 150 m is assumed, covering an area of 10 km^. 

The angular resolution of individual events is found to be about 0.3° near the energy threshold, improving to below 0.1° at higher 
energies. The relative energy resolution is 20% at the threshold and improves to 10% at higher energies. Several parameters for 
gamma hadron separation are described. With a combination of these parameters, 80% to 90% of the hadronic background can be 
suppressed, while about 60% of the gamma-ray events are retained. The point source sensitivity to gamma-ray sources is estimated. 
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O optimistic assumptions, and a 100 km"^ array, a sensitivity of about 1 x lO"'"* erg s"' cm""^ can be achieved (at 100 TeV) 
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using conservative assumptions, to be about 8x10 '-^ erg s ' cm ^ at 100 TeV gamma-ray energy for a 10 km^ array. With more 
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■ Even in the former case the detector is sensitive enough to measure the continuation of currently known gamma-ray spectra 
into the ultra high energy domain. Due to its large field of view of 0.6 sr it also offers a great potential for the discovery of new 
gamma-ray sources at the so far largely unexplored energies of 100 TeV and above. 

fCeywords: 
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Reparation 



1. Introduction 

Ground-based gamma-ray astronomy has become a success- 
ful and well-established branch of modem astronomy and of- 
fers unique possibilities to study the extreme conditions in the 
non-thermal universe. More than hundred sources of very high 
energy (VHE, 30 Ge V to 30 TeV) gamma-ray emission are cur- 
rently known, and detailed spectral, morphological and tempo- 
ral studies provide a wealth of information on the underlying 
acceleration and radiation mechanisms of these sources (see 
e.g. Ctl, lH or H for reviews). 

However, at even higher energies, in the ultra high energy 
band (UHE, Ey > 30 TeV), not many sources could be discov- 
ered so far. The main challenge at these energies is the very low 
photon flux emitted even by strong gamma-ray sources. Cur- 
rent instruments for gamma-ray observations, like the 3'^'' gen- 
eration Cherenkov telescope arrays [4i, are usually optimised 
for a good sensitivity in the VHE band. While their effective 



areas of around 10^ m^ are well suited for this energy band, 
they are too small for sensitive observations in the UHE band: 
For example, in 23 hours of observations of the Crab Nebula 
with H.E.S.S., only four gamma-ray events with energies above 
30 TeV could be found ||5|]. 

On the other hand, observations at ultra high energies can 
help to answer important open questions in high energy as- 
tronomy. One of the most important scientific goals of UHE 
gamma-ray astronomy is the search for the sources of charged 
cosmic rays at PeV energies (PeVatrons) It is believed 
that cosmic rays at least up to the so-called knee of their en- 
ergy spectrum at around 4 PeV are of Galactic origin but 
so far no conclusive evidence could be found for sources or 
source classes capable of accelerating particles to that energy. A 
clear evidence could be a gamma-ray signal above 100 TeV fl. 
However, no such signal could be observed so far, possibly due 
to the lack of sensitive observation^ The gamma-ray spectrum 
at ultra high energies can also be used to distinguish effectively 
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' In fact, also a high-energy neutrino signal could be an unambiguous indi- 
cation for cosmic ray acceleration, but so far none has been found either j^. 
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between leptonic and hadronic cosmic ray accelerators - a noto- 
riously difficult task with VHE band observations alone - since 
the cross section for leptonic gamma-ray production decreases 
strongly in the UHE regime (Klein-Nishina regime). Further 
motivations for UHE gamma-ray astronomy can be found in 

Jilii. 

To overcome the challenge of small event numbers in UHE 
observations, very large effective areas (in the order of many 
km^, ideally 100 km^ or more) and long exposure times per 
source are required. A powerful concept that can achieve 
these requirements are arrays of non-imaging, wide-angle air 
Cherenkov detectors. In these detectors each station samples 
a part of the Cherenkov light front without any image infor- 
mation. Typically, the detector stations contain only one or a 
few large photomultiplier tubes (PMTs) directed towards zenith 
with a light concentrator on top to increase the light sensitive 
area and to restrict the influence of stray Ught from large zenith 
angles. The modest light collection areas of (9(1 m^) (com- 
pared to (9(100m^) in modern imaging Cherenkov telescopes) 
are sufficient since at ultra high energies each incident gamma 
photon produces a very large air shower and thereby a rela- 
tively strong Cherenkov light pulse. The inexpensive design of 
the individual stations and large inter-station distances make the 
realisation of large effective areas feasible. Additionally, indi- 
vidual sources receive on average a higher exposure than from 
Cherenkov telescopes, since a large part of the sky is monitored 
simultaneously (up to 1 sr). 

Experiments of this kind have first been developed for studies 
of charged cosmic rays, e.g. at Yakutsk ifTlll or at the Buck- 
land park observatory in Australia [T?]. More recent detec- 
tors employing this method are the BLANKA 11311 and the 
TUNKA 11411 experiments, the latter still being operational, 
again with the focus on charged cosmic rays. The only non- 
imaging Cherenkov detector array specifically designed for 
gamma-ray astronomy has been the AIROBICC detector built 
by the HEGRA collaboration [15]. With an instrumented area 
of only 3 x lO'^m^ this array was however too small to detect 
gamma-ray sources above its energy threshold of about 15 TeV. 

Currently, the HiSCORE detector is being designed as a new, 
powerful instrument for gamma-ray astronomy at ultra high en- 
ergies. With an effective area of 10 to 100 km^, fast read-out 
electronics and advanced event reconstruction techniques it is 
expected to achieve a gamma-ray sensitivity that allows a de- 
tailed study of the UHE part of the spectrum of the currently 
known sources and the potential detection of many new sources 
in the UHE regime. By its design it will also be a powerful 
instrument for the measurement of charged cosmic rays, with 
a focus on the very interesting energy range around the knee. 
Details of the detector can be found in UL^. 

To achieve the large effective area at a reasonable effort, an 
inter-station spacing of 100 m to 200 m is foreseen (to be com- 
pared with, e.g., approximately 15 m for the individual stations 
that made up the AIROBICC array). The resulting low number 
of data channels per event requires advanced air shower recon- 
struction techniques that can work with a sparse sampling of 
the Cherenkov Ught front. In this paper, the event reconstruc- 
tion developed so far for this detector will be presented. All 



reconstruction techniques are tested using Monte Carlo simula- 
tions and the results are used to calculate the expected accuracy 
of the reconstruction. 

In section|2]the air shower and detector simulations are intro- 
duced along with the extraction of signal parameters from the 
data. Section [3] outlines the reconstruction of shower core po- 
sition, direction and energy of the primary particle and the ver- 
tical position of the Cherenkov light maximum (shower depth). 
The particle separation algorithm is presented in section |4] The 
results are used in section |5] to calculate the gamma-ray flux 
sensitivity of the instrument. 

2. Detector simulation and measured quantities 

2.1. Air shower simulations 

The simulation chain used for this work consists of the air 
shower simulation cod^ CORSIKA [17J, the detector simu- 
lation sim_score ifisll and the event reconstruction reco^core. 
About 145,000 air showers with energies between 10 TeV and 
5 PeV have been generated following a dN/dE oc spectrum. 
Primary particle species used include gammas, protons and he- 
lium with about 25,000 events each and nitrogen and iron with 
about 35,000 events each. Event core locations are scattered 
uniformly over the array and up to 350 m beyond the detector 
perimeter Particle directions are randomised with a maximum 
zenith angle of 30°. The CORSIKA lACT package [19] has 
been used to obtain the Cherenkov light pulses at the detector 
level. All simulations are carried out using the US standard at- 
mosphere [20]. The impact of varying atmospheric conditions 
on the results has not yet been studied. 

2.2. Detector simulation 

The current simulations of HiSCORE assume a regular 
square grid of 22 x 22 detector stations with an inter-station 
distance of 150 m at sea level, resulting in an instrumented area 
of 9.922 km^. (The effects of a different station spacing and al- 
titude are discussed in [23[, and further investigations into an 
optimised layout, possibly using varying stations spacing, are 
on-going.) 

Each detector station contains four 8" PMTs equipped with 
Winston cone light concentrators, and has a total light sensitive 
area of 0.5 m^. The signals of all four channels are added up and 
recorded by a 1 GHz analogue-digital converter Each station is 
triggered independently by signals greater than about 180 pho- 
toelectrons. Only triggered stations and stations that have at 
least one triggered neighbour station are read out. The sam- 
pled signals and the corresponding time stamps are transmitted 
to the central data acquisition system digitally, using standard 
technology like ethernet or a mobile phone network. 

The detector simulation sim^core calculates the signals seen 
by each station, taking into account atmospheric absorption, 
angle- and wavelength dependent transmission of the Winston 



^Using CORSIKA version 6.735 witli tlie QGSJET module version Olc.f 
(2lll for high energy hadronic interactions, the GHEISHA module 1. 2^1 for 
lower energies, and EGS4 for the electromagnetic component. 
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Figure 1: Typical detector signal with illustration of the signal parameters in- 
troduced in section 12.31 The signal is extracted from a simulated 509 TeV 
gamma-ray event, from a detector station about 200 m away from the shower 
core (summed signal of all four PMT channels). The average noise level of 
90p.e./ns has been subtracted. The x-axis shows time relative to an arbitrary 
reference. 



cones, the PMT quantum efficiency, PMT afterpulses, the signal 
shaping by the system's response function, and the noise in the 
system, which is dominated by light from night sky brightness 

\m 



2.3. Data processing 

The signals are preprocessed for the reconstruction frame- 
work by extracting the following parameters (see also figure[TJ: 

• The time of the signal peak (peak time), defined as the time 
bin with the highest entry 

• The time where the signal reaches 50% of the peak (edge 
time), linearly interpolated between adjacent time bins 

• The rise time, defined as the time in which the signal rises 
from 20% to 80% of its maximum value (interpolated) 

• The full width at half maximum (FWHM) of the signal 
(signal width) 

• Signal size, estimated by integrating the signal from 5 
nanoseconds before the peak time to 15 nanoseconds af- 
ter the peak time 

The full sampled waveforms of the signals are stored addi- 
tionally for later use in the reconstruction. 

2.4. Quality cuts 

To select a sample of good events, basic quality cuts are ap- 
plied to the data (acceptance cuts). In this study, events are used 
if they trigger at least three stations, have a reconstructed core 
position within the array and a reconstructed zenith angle of no 
more than 25°. This last requirement determines the effective 
field of view to 0.6 sr. Loosening these cuts would deteriorate 



the reconstruction accuracy, but also improve the event statis- 
tics, and may be beneficial in some cases - this is however not 
studied here. On top of the acceptance cuts, gamma-ray cuts 
are used to select a sample of gamma-ray like events. These 
cuts are described in section |4] 

3. Event reconstruction 

The event reconstruction starts with the reconstruction of the 
shower core position and the direction of the primary particle 
as pure geometric parameters. Subsequently, the energy and the 
shower depth are reconstructed using lookup tables. The recon- 
structed values are then used for the gamma hadron separation 
(see sectionlHi. 

In the following, the resolution of each reconstructed quan- 
tity will be given as the value at which 68% of the events are 
contained. All gamma-ray events from the dataset described in 
the previous section that survive the acceptance and the gamma- 
ray cuts are used to calculate the resolution. Generally, all res- 
olutions improve with the number of triggered stations, and 
hence with energy of the primary particle, and will therefore 
be given as function of energy. Only data points above 50 TeV 
are presented, since the reconstruction is not yet optimised for 
lower energies and deteriorates quickly in that regime. The 
question of the energy threshold will be discussed briefly in sec- 
tion lO 

3.1. Shower core position 

The most straightforward method to reconstruct the position 
of the shower core (the intersection of the shower axis with the 
detector level) is a centre of gravity calculation using signal 
sizes of all triggered stations, as used e.g. in the AIROBICC 
experiment 11511 . This is a robust method that usually gives a 
good approximation of the core position even if only a few sta- 
tions have triggered. 

If signals from at least five stations are available, a better 
core position can be obtained by fitting the lateral light density 
function (LDF) to the recorded intensities, as suggested in |25]. 
The LDF is being parameterised as an exponential function near 
the shower core and a power law at larger distances, with the 
break between the two at c^df ~ 120 m 11311 : 



LDF(r) = 



P exp(d r) for r < cldf 
Qr'' for r > cldf 



with 
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r(x,y) 

P exp(dcLDF) 
(cldf)'' 



+ y'^ 



(1) 

(2) 
(3) 



The free parameters of the fit are the absolute normalisation 
P, the inverse decay-length of the exponential function d, the 
power law index k and the position of the shower core (x,y). 
If at least six detector signals are available, ciof can be a free 
fit parameter as well, which has been found to slightly improve 
the core resolution. 
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Figure 2: The LDF of two 750 TeV gamma-ray events with different shower 
depths. The distributions are fitted with equation \l\ The dashed lines denote 
the exponential pait of the LDF, the solid lines the power law part. 



The data is fitted using an iterative ;^f^-minimisation. The 
data point with the largest deviation from the fit (in terms of 
cr) is removed, and a new fit is conducted using the reduced 
data sample. This procedure is repeated until there are no more 
outliers further away than 1 a from the fit, or the number of 
remaining data points drops below five. The fit is considered 
successful if it converges numerically, the resulting value is 
not used for further evaluation. The same procedure is used for 
all following fits as well. 

After the reconstruction of the shower direction (see next sec- 
tion), the station positions can be transformed into the shower 
plane (the plane perpendicular to the shower axis), and the core 
position fit can be repeated using the new coordinates, which 
yields a slight improvement especially for events with large 
zenith angles. Figure |2] shows two examples of the LDF fit- 
ted to simulated data (two-dimensional projection). The use of 
the LDF for the energy and shower depth reconstruction is dis- 
cussed in sections [331 and f3A\ 

Figure |3] presents the resolution of the shower core position 
reconstruction for the centre of gravity method and the LDF fit 
method. The resolution is about 35 m near the threshold and 
improves to less than 10 m at higher energies. 

While the centre of gravity method can (per definition) only 
reconstruct core positions within the array, the LDF fit can also 
be used to reconstruct core positions outside of the array. How- 
ever, this method gets inaccurate quickly when the shower core 
moves away from the array border At 1 PeV, the resolution 
deteriorates from 5 m for contained events to 15 m for core po- 
sitions up to 150 m away from the detector perimeter, and to 
about 35 m for showers between 150 m and 300 m away. 

A better method for these showers is the usage of the signal 
width, which can be used up to far distances from the array as 
long as enough light is received in the detector to accurately 
measure the signal widths. This method has successfully been 
tested by the AIROBICC collaboration ll26ll and is also used at 
the TUNKA detector ll4ll to increase its effective area. This 
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Figure 3 
struction 



Resolution (68% containment) of the shower core position recon- 
with the two methods described in section IsTI 



method has however not been pursued in this study, and only 
events that are reconstructed to be inside of the array are used 
for further analysis. 

3.2. Direction of origin 

A coarse estimate of the shower direction can be obtained 
by fitting a plane to the measured light arrival times0. If signals 
from only three stations are available, only this estimate is used. 

If at least four signals are available, the direction reconstruc- 
tion can be improved significantly by taking into account the 
curvature of the Cherenkov light jfront. For this, the light ar- 
rival time model developed by ll27ll for timing stereoscopy at 
Cherenkov telescopes has been adapted for an array detector. 
The expected time delay toer at a given detector - with respect 
to the arrival time at the core position, fo - is parameterised 
as function of the detector position, the height of the shower 
maximum z, and the direction of the shower axis (given by the 
zenith angle and the azimuth angle <p). The detector position 
is given relative to the shower core in polar coordinates, using 
its distance r and the azimuth angle ipoei- 



tOet - toetir, 4>Det, Z, 9, (fl) + tQ 



(4) 



Using r and (poet from the previous reconstruction of the 
shower core position, the function can be fitted to the measured 
arrival times, which yields the height z and the shower direc- 
tion. 

The functional dependence is derived by integrating the light 
path from the point of Cherenkov light emission to the detector, 
using 

Tj{h) = 1+770 exp{-h/ho) (5) 

for the height dependent refractive index of air (ho = 8 km, 
770 = 2.76 X 10""^). The integration yields 



^^Either peak times or edge times can be used. Using the edge times results 
in a slightly better accuracy, as they can be measured with higher resolution. 
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Figure 4: Resolution (68% containment) of tlie direction reconstruction. Tliree 
different scenarios are assumed for the time synchronisation of the signals be- 
tween the stations: An ideal synchronisation, a synchronisation with a gaussian 
en'or with cr = 1 ns and one with cr = 2 ns. 



tdet(k,Z) 



+ to 



(6) 



with 



k = k(r,<p',z,e) ^r^ + 



- 2rztan(0) cos((p') (7) 



with (/)' = (poet - 4> (see H28ll for details on the derivation). 

This model is fitted to the measured arrival times (edge times) 
with the height z and the direction described by 6 and (p as free 
parameters. The fit results in an accurate reconstruction of the 
particle direction and an estimate for the height of the shower 
maximum. Other methods to reconstruct the shower maximum 
are described in section l34l 

A good angular resolution can only be achieved if the signals 
of the different stations can be synchronised to each other with 
high precision, which is a challenging experimental task given 
the large distances between the stations. To examine the im- 
pact of a non-ideal time synchronisation, the simulated signals 
have been shifted randomly in time using a Gaussian distribu- 
tion with a width of <j. Figure |4] shows the angular resolution 
as function of energy for a perfect synchronisation and for jit- 
ters of cr = 1 ns and cr = 2 ns. The angular resolution is given 
as the 68% containment region of the angular distance between 
reconstructed and true direction of individual events. 

The results show that the accuracy of the direction recon- 
struction is limited by the time synchronisation even if the jitter 
is as low as 1 ns. Therefore, a time synchronisation accuracy of 
at least 1 ns, preferably better, should be aimed at in the detec- 
tor development. In the following, a time synchronisation with 
a jitter of cr = 1 ns will be assumed. With this assumption the 
resolution of the direction reconstruction ranges from 0.4° near 
the threshold down to 0. 1° at higher energies. 
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Figure 5: The resolution (68% containment) of the energy reconstruction as 
described in section [53] 



3.3. Energy 

In general, the amount of Cherenkov light generated is pro- 
portional to the energy of the primary particle. Since only a 
small fraction (about 0.002%) of the total Cherenkov light ar- 
riving at the observation level is captured by the HiSCORE de- 
tector stations (due to their large spacing and their small light 
collection areas), the fitted LDF is used for an estimate of the 
total Cherenkov light. The actual energy reconstruction is com- 
plicated by the fact that the light distribution on the ground de- 
pends strongly on the shower depth: The lower a shower max- 
imum occurs in the atmosphere, the steeper the LDF will be, 
with more light closer to the shower core and less further away 
(see also figure|2]i. 

The simulations carried out show empirically that the impact 
of the shower depth on the light distribution at detector level is 
minimal at a distance of around 220 m from the shower core. 
Therefore, the fitted light intensity at 220 m, i.e. LDF(220m), 
is used for energy reconstruction. Simulated gamma-ray events 
are used to generate a lookup table of log(LDF(220m)) to 
log(£') with six bins per decade. The energy is reconstructed 
from this lookup table using interpolation. Figure |5] shows the 
resulting relative energy resolution. It improves from about 
30% near the threshold to about 10% at higher energies. The 
selection bias that usually leads to an overestimation of particle 
energies near the threshold is negligible above 50 TeV. 

The energy reconstruction is optimised for gamma-ray 
events. Hadronic cosmic rays produce less Cherenkov light 
than gamma-rays with the same energy, and heavier hadronic 
primaries produce less light than lighter ones. The energy of 
cosmic rays, especially heavier nuclei, is therefore always un- 
derestimated. To derive a cosmic ray energy spectrum, an as- 
sumption about the composition must be made. 

3.4. Shower depth 

The shower depth is defined as the atmospheric depth (mea- 
sured from the top of the atmosphere along the shower axis) of 



5 



05 



c 

(/) 



45 
40 
35 
30 
25 
20 
15 
10 
5 

-5 

50 100 150 200 250 300 350 400 
Core distance [m] 

Figure 6: The signal widths of the same two gamma-ray events as in figure|2] 
after summing up signals from stations within 50 m bins. Linear fits are applied 
in the core distance range from 150 m to 400 m. The errorbars are estimated on 
basis of the signal size after stacking. 



d 

" d 


epth = 
epth = 


= 495 
= 718 


2 

g/cm 

g/cm 






















































































































J 



















































E 



00 

c 
o 

o 

03 
03 



03 
Q 



90 
80 
70 



30 



1 1 

X 


1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 

Timing + 
LDF X 
widths O 
combined • 


o 
• 


o 






• X S o 

• X * 






X o 
• X ^ 

• • 


® O ^ 


.... 


1 1 1 1 1 ' ' ' 


• + • 
• 



10" 



lO'' 
Energy [TeV] 



10" 



Figure 7: Shower depth resolution of proton events for the three methods de- 
scribed in section lT4l and the combination of the three methods, versus the true 
(simulated) energy. 



the maximum of Cherenkov light emission, which (almost) co- 
incides with the maximum of relativistic particles in the shower. 
Although not of immediate interest, it is an important quantity 
for the gamma hadron separation (see section |4) and the deter- 
mination of the mass (in case of a hadronic primary). Three 
different (partly correlated) methods are used here for the esti - 
mation of the shower depth: 

Arrival time method. As shown in section 13.21 the height of 
the light emission has an impact on the shape of the arrival 
time distribution over the array. Generally, the delay of the 
photon arrival at stations at a given distance from the shower 
core increases with shower depth. Using the peak times of the 
Cherenkov light signals, the height can be inferred from the fit 
to the function toet (eq- ISll. In practice, the reconstruction has 
proven to be more accurate if the edge times are used instead 
of the peak times. A lookup table (fitted value of z to shower 
depth) is used to correct for the shift introduced by using the 
edge times. 

LDF slope method. As illustrated in figure |2] the Cherenkov 
light of showers with a large shower depth is concentrated 
closer around the shower core position. Usual parameters for 
the measurement of this feature are the slope of the inner part 
of the LDF or the ratio of the LDF values at a small and a large 
core distance S IM El H Q . Since the slope of the in- 
ner part of the LDF is not well defined in HiSCORE due to 
the large station spacing, the ratio LDF(50m)/LDF(220m) is 
used. Tests with other ratios, e.g. LDF(50m)/LDF(150m), 
have yielded similar results (the best ratio also depends slightly 
on the energy range). Lookup tables are used to derive the depth 
from the measured ratio. 

Signal widths method. The width (or duration) of the 
Cherenkov light pulse - defined here as full width at half max- 
imum, FWHM - increases with core distance r. Different pa- 



rameterisations for width(r) have been suggested, e.g. a power 
law (y^ or an exponential function |25|, but for the current 

H) has 



study a simple linear function (as also used by |134 
been found to describe the simulated data best. Close to the 
shower core the pulse is shorter than the time resolution of the 
simulated PMT, so that this functional dependence is only visi- 
ble in detectors further than about 150 m away from the shower 
core. 

At a given core distance, the signal width increases with 
shower depth, as observed already experimentally by ||34 35 ] 
and studied in simulations by (36ll . It can therefore be used 
for the determination of the shower depth, as done e.g. in the 
TUNKA analysis |37]. The sensitivity to the shower depth in- 
creases with core distance, but the determination of the width 
becomes more and more inaccurate at large core distances due 
to the decreasing signal sizes, so a trade-off must be found. 

To improve the measurement of the signal FWHM at large 
core distances, all signals within a certain distance interval from 
the shower core are summed up before determining the sig- 
nal width (signal stacking). As the number of detector stations 
available in a distance interval is proportional to the distance to 
the shower core, this procedure makes it possible to sample the 
width up to high distances from the shower core. The distance 
intervals used here are 50 m wide. Figure |6]shows the resulting 
widths distribution of two events with different shower depths, 
along with linear fits in the core distance range from 150 m to 
400 m. 

The depth reconstruction is done here using the fitted FWHM 
at 300 m in combination with lookup tables. 



Combination of methods and performance. To combine the 
three methods, the average of all methods that return a success- 
ful depth reconstruction is used (a method may be unsuccessful 
if the required parameter distribution cannot be fitted due to 
missing or contradictory data). Figure Q shows the resolution 
of the three individual methods and the average for a sample of 
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proton events. Since all methods require a good determination 
of parameter distributions up to large core distances, a depth 
estimate is difficult near the threshold where only a few stations 
produce a usable signal. Towards PeV energies, the resolution 
improves to below 40 g/cm^. 

The simulations show that a particle dependent bias exists in 
all three depth reconstruction methods. Therefore, the depth 
reconstruction can only be optimised to one particle species; 
for the current study, it is calibrated to protons. The causes 
for this behaviour and the resulting potential for gamma hadron 
separation are discussed in section |4T| 

4. Gamma hadron separation 

The objective of gamma hadron separation is the statistical 
separation of gamma-ray induced air shower events among the 
(more or less) constant and isotropic abundance of (charged) 
cosmic ray eventfl 

Since the observable differences between photonic and 
hadronic air showers are rather subtle, it is traditionally a diffi- 
cult task to distinguish them with any air Cherenkov detector, 
but particularly with non-imaging detectors. Nevertheless the 
separation is possible to some extent, as will be shown. 

It should be noted that small gamma-ray emission regions, 
especially point-like sources, can be detected without any 
gamma hadron separation, as they emerge from the isotropic 
flux of background as a localised excess. Nevertheless, a good 
gamma hadron separation improves the significance of a detec- 
tion, or in other words, enables the detection of weaker gamma- 
ray sources during the same observation time. 

In the following, different methods for particle separation 
will be discussed qualitatively. The combination of separation 
parameters and the resulting quantitative performance will be 
discussed section |4~: 



4.1. Separation using depth reconstruction bias 

The simulations show that all three described depth recon- 
struction methods exhibit a particle dependent bias. With the 
LDF and the widths method, the depth of showers induced by 
heavier hadrons (e.g. iron) is systematically overestimated by 
60 to 100 g/cm-, while the depth of photonic showers is under- 
estimated by 20 to 40g/cm^. Since the lookup tables are done 
with a proton sample, no bias exists for protons. If the arrival 
time method is used, the offset between gamma-rays and iron 
nuclei, the most extreme cases among the simulated species, 
amounts to only 20 to 50 g/cm^. 

The reason for these offsets are differences in the shower de- 
velopment below the shower maximum that depend on the type 
of primary particle. Air showers induced by hadrons, especially 
heavier nuclei, contain large numbers of secondary hadrons. 
These particles do not loose energy as rapidly as particles from 



"^Since the flux of electrons decreases more rapidly with energy than the one 
of nuclei, electron-induced air showers are not a relevant cause of background 
for UHE gamma-ray astronomy. Therefore, only air showers of hadronic origin 
must be filtered out. 



the electromagnetic cascade and therefore penetrate deeper into 
the atmosphere. Even at low altitudes, they constantly re- 
fuel (e.g. by pion decay) the electromagnetic component of the 
shower, which results in more Cherenkov light closer to the 
ground. Since the LDF and widths methods are based upon the 
full sampled Cherenkov light pulses, they are sensitive to all 
stages of the shower development. The Cherenkov light from 
low altitudes shifts the reconstructed shower maximum to there, 
i.e. to larger depths. The arrival time method, on the other hand, 
is only sensitive to the Cherenkov light emitted at a specific 
point of the shower development (e.g. the shower maximum, if 
the peak times are used) and is therefore not influenced by the 
differences between the shower types. 

There are several ways to exploit these facts for particle sep- 
aration. Here, the difference between the depth reconstructed 
with the widths method, X^'idih, and the arrival time method, 
Xiiming, is used. The centre of the distiibution of X^i^ih -X,i,„i„g is 
around zero for protons (since the lookup tables for both meth- 
ods are done with a proton sample), below zero for gammas and 
above zero for heavier nuclei. Although the distributions over- 
lap significantly, a particle separation can be achieved to some 
extent. 

It should be noted that the simulations indicate that the offsets 
between the reconstructed depth values are most pronounced 
for a detector at sea level and almost vanish at an altitude of 
2000 m, probably because the later stages of the shower, which 
are sensitive to the hadron content, do not even develop above 
this altitude. For the same reason, events with higher ener- 
gies (starting at low PeV energies) show smaller offsets, as 
their shower depth is larger on average and the later stages 
of the shower can no longer develop above the detector level. 
This may explain why such a bias has not yet been noted by 
other non-imaging air Cherenkov arrays (e.g. AIROBICC 113 ill . 
BLANCA 1 13] or TUNKA |37]): Either their altitude was too 
high, or the focus was on too large energies, or both. 

4.2. Separation using shower depth versus energy 

The depth of an air shower is largely determined by the alti- 
tude of the first interaction of the primary particle in the atmo- 
sphere. The altitude of the first interaction, in turn, depends on 
the cross section of the reaction between the primary particle 
and air molecules. The shower depth increases logarithmically 
with energy and is larger for photonic than for hadronic events. 
Additionally, heavier nuclei result in a smaller shower depth 
than lighter nuclei at the same energy. Therefore, the shower 
depth (in combination with the previously reconstructed parti- 
cle energy) can be used as an indication of the particle type. 
This method is widely used in the field of air shower detectors, 
usually to derive the mean mass of charged cosmic ray particles 
(see e.g. 101). 

The particle identification is complicated by the fact that the 
depth of hadrons, especially heavier hadrons, is usually overes- 
timated by the described depth reconstruction methods. To keep 
the bias on the reconstructed shower depth as small as possible, 
only the arrival time method is used here. Additionally, the 
energy of hadronic particles is underestimated. Both offsets de- 
crease the observable differences between the particles. The 
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Figure 8: Median reconstructed depth (arrival time metliod) versus recon- 
structed energy for gammas, protons and iron nuclei. The errorbars denote 
statistical variation (68% containment), i.e. include both the shower fluctua- 
tions and the uncertainty of the reconstruction. Values for protons are shown at 
10%, for iron nuclei at 20% higher energies for clarity. 



remaining subtle differences between nucleonic and gamma- 
induced air showers are washed out to some extent by the in- 
trinsic fluctuations of the air shower development. Neverthe- 
less, some separation of gammas and hadrons can be achieved. 

Figure |8] shows the median reconstructed depth values for 
different particles as function of reconstructed energy. The me- 
dian reconstructed depth for photonic events is parameterised as 
function of energy. The deviation of a measured value from this 
expected depth is used as particle separation parameter While 
the overlap between the distributions for gammas and protons 
is considerable, heavier nuclei can be identified and discarded 
rather efficiently. 

4.3. Separation using signal rise time 



As noted in section 14.11 air showers of hadronic origin con- 
tain more secondary hadrons which refuel the electromagnetic 
cascade down to low altitudes. Since these secondary hadrons 
move faster than the speed of light in air, the Cherenkov light 
produced by the corresponding electromagnetic cascade ap- 
pears in the detector before the bulk of Cherenkov light gen- 
erated from particles near the shower maximum. Therefore, 
hadronic events can be identified by some "early light" before 
the peak of the signal. 

Here, this feature is detected by using the signal rise time (see 
section 12.31 1. For hadrons (especially heavier ones), a longer 
rise time is expected. The effect is most pronounced near the 
shower core, since at larger core distances the light from low 
altitudes must cross a significantly longer distance through air 
and no longer appears before the main peak. Therefore, only 
the signal from the "central station" (the station closest to the 
shower core) is used. Figure |9] shows the distribution of this 
parameter for various particles and the main energy range of 
HiSCORE. While again the distributions overlap considerably, 
a separation is evident. 
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Figure 9: Distribution of the signal rise time of the central detector station for 
gammas, protons and iron nuclei with (true) energies between 100 TeV and 
IPeV. 



The rise time has been suggested previously as particle sep- 
aration parameter for Cherenkov telescope data (see e.g. [13 811 ). 
but due to its close correlation with the image shape recorded in 
the telescopes it was considered not useful |39i]. It seems how- 
ever that for non-imaging detectors the rise time can be a useful 
parameter for gamma hadron separation. 



4.4. Separation using Cherenkov light spectrum 

The Cherenkov light emitted by hadronic air showers at low 
altitudes can not only be identified by its time signature, but 
also by the observable light spectrum. Since UV light (around 
250 nm) is strongly absorbed by air, its presence at the detector 
level is an indication of light emitted at low altitudes, and thus 
of hadronic events. 

The use of this method has been suggested previously by 
ll40ll . and applied to simulations of imaging Cherenkov tele- 
scopes by i4lll and i38fl . Although a signature of the particle 
type could be found in the amount of measured UV light, these 
studies concluded that other, easier to measure, parameters con- 
tain the same or even better information. 

However, these studies concentrated on imaging Cherenkov 
telescopes, and it seems plausible to assume that the spectral 
information may be more useful in a non-imaging detector (in 
a similar way as the rise time, see section l431 l. Therefore, the 
effect of a spectral measurement for HiSCORE was tested by 
adding an additional UV sensitive channel to the standard de- 
tector stations in the simulation. The amount of detected UV 
light is normalised to the light detected in the standard detec- 
tion channel, yielding the UV light ratio. 

The simulations show that almost no particle separation can 
be achieved by this method, as two contrary effects cancel out 
each other: On the one hand, the UV light ratio contains a sig- 
nature of the shower depth. Since air showers of photonic origin 
penetrate on average deeper into the atmosphere than hadronic 
showers (see e.g. figure |8]l, they contain more particles close 
to the detector level, which increase the UV light ratio. On 
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the other hand, hadronic events exhibit more deeply penetrat- 
ing secondary hadrons which produce some early light close to 
the detector (see explanation in sections|4T|and|43]l and by that 
a surplus of UV light. 

A correction for the (reconstructed) shower depth can be 
used to isolate the latter effect, but introduces additional uncer- 
tainties. The remaining differences between different particle 
species are much smaller than in the previously discussed meth- 
ods. Therefore, the additional effort of a spectral measurement 
does not seem to be justified. 

4.5. Combination of methods and performance 

To rate the performance of the gamma hadron cuts, the qual- 
ity factor QF is used: 



2.5 



(8) 



where ey and ebg are the survival probabilities for gamma-ray 
and cosmic ray (background) events, respectively. In the (com- 
mon) case of background-dominated observations, the quality 
factor has a linear impact on the instrument sensitivity. 

The survival probabilities for each particle species can be de- 
termined by applying the gamma hadron cuts to the simulated 
events. However, since all described particle separation meth- 
ods discriminate heavier hadrons more efficiently than lighter 
ones, the correct calculation of e^g must assume a realistic com- 
position of the background events. For this, the polygonato 
model by 14211 is used, which yields parametrisations of the en- 
ergy spectra of all elements up to Z = 92. The survival proba- 
bility for events of each particle species are calculated using the 
particle separation power of the most similar simulated element 
(iron, nitrogen, helium, hydrogen). 

Additionally, the non-ideal energy reconstruction must be 
taken into account: Since the energy of cosmic rays is underes- 
timated on average (see also section 13.3b . the background at a 
given reconstructed energy comprises of cosmic ray events with 
a higher true energy. In order to compare and etg, both are 
always given for reconstructed energies. 

Several strategies can be used to combine the three particle 
separation methods introduced in sections l4n to 14. 3l (the spec- 
tral method discussed in section l4~4l is not used). Here, the cuts 
on each parameter are applied successively to each event. The 
(energy dependent) cut values for each parameter are adjusted 
to leave about 80% of the gamma-ray events at each stage. The 
resulting total gamma-ray survival probability ranges from 
50% to 60%, while ei,g improves from 30% at lower energies 
to less than 10% at higher energies (see figure [TOt. The cor- 
responding quality factor improves with energy from 0.9 near 
the energy threshold to about 2.0 at PeV energies. A more so- 
phisticated way to combine the three cuts, e.g. a multivariate 
analysis, may improve the performance slightly (see e.g. [43] 
for a corresponding work for imaging Cherenkov telescopes). 

It should be noted that, depending on flux level, observation 
time and other factors, there may be an energy above which the 
sensitivity is limited by the number of signal events rather than 
by the background (background-free regime). Obviously, the 
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Figure 10: Gamma-ray and cosmic ray (background) survival probabilities, and 
quality factor, after the application of the described combination of gamma 
hadron cuts, versus the reconstructed energy. The uncertainties of the survival 
probabilities are about the size of the markers or smaller 



gamma hadron separation should be adjusted at these energies 
in order to retain as many gamma-ray events as possible. 

Near the energy threshold, up to a few hundred TeV, al- 
most no gamma hadron separation can be achieved with the de- 
scribed strategy. The reason for this is the poor accuracy of the 
depth reconstruction in this regime, which deteriorates the sep- 
aration of the methods described in sections WA\ and W2\ How- 
ever, using the rise time method (section 14.3b alone, a quality 
factor of 1.2 to 1.3 can be achieved from 50 TeV on. At higher 
energies, the rise time alone is less effective than the described 
combination of parameters. In practice, the most appropriate 
gamma hadron separation scheme should be selected after the 
reconstruction of the energy. 

5. Detector sensitivity 

5.7. Effective areas 

The effective areas of the detector are given by the instru- 
mented area A (10 km^ for these simulations) and the average 
probabilities for events to survive both the acceptance (eacc) and 
the gamma hadron cuts (e^ or ebg, depending on the type of par- 
ticle): 



^eff,y/hg 



£y/bg 



(9) 



Like the survival probabilities eacc, £y and ebg, they are a func- 
tion of energy and particle. Figure [TT] shows the effective areas 
as determined from the simulations, for all five simulated parti- 
cle species. As expected, the effective areas for heavier nuclei 
are reduced efficiently by the gamma hadron separation. 

From the effective areas it can be seen that gamma-rays 
can be detected from energies of about 30 TeV. At about 
50 TeV, the effective area for gamma-ray detection reaches 
3 km^, roughly 50% of its maximum value. 

For the study of charged cosmic rays, the analysis will be 
performed without the gamma hadron separation. In this case. 
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Figure 1 1 : Effective areas for all five simulated particle species after gamma 
hadron separation. Errorbars denote statistical uncertainties estimated from the 
expected Poissonian variance in each bin. 



the effective areas are equal to the instrumented area (lOkm^ 
in the simulated case) above a particle dependent energy (about 
120 TeV for protons, about 300 TeV for iron nuclei). 

5.2. Point source sensitivity 

The point source sensitivity of the detector is one of the 
key performance figures and can be compared with the flux of 
gamma-ray sources. It is estimated here as the minimal flux 
of gamma-rays that is needed to detect a point-like gamma-ray 
source with 5cr above the uniform background of cosmic ray 
events and at least 50 gamma-ray events (detection criteria). 

The rate of background events is calculated using the flux 
parameterisations from [42] together with the eff'ective areas 
shown in figure [TT] The systematic underestimation of the cos- 
mic ray energy is taken into account and reduces the flux seen 
at a given reconstructed energy by about a factor of two. 

To calculate the number of gamma-ray events required for a 
detection, the rate of background events in the source region 
must be known. For the estimation of the sensitivity, a circu- 
lar source region with a radius equal to the angular resolution 
shown in figure |4] is assumecjfl For a conservative estimate, 
an inter-station time synchronisation of 1 ns is assumed. For a 
more optimistic estimate, an ideal time synchronisation is as- 
sumed. Since the sensitivity is proportional to the square root 
of background events, it depends linearly on the angular reso- 
lution. 

Figure [12] shows the calculated rates of cosmic ray events for 
the whole field of view and within the source region. The values 
for the eff'ective areas and the angular resolution are inter- and 
extrapolated where needed. 

To determine the background level that needs to be subtracted 
from the number of events in the source region, several methods 



Figure 12: The integral cosmic ray event rates seen by HiSCORE before and 
after gamma hadron separation. The third line demonstrates the effect of the 
underestimated reconstructed energy for hadronic events. Lines 1 to 3 show 
the cosmic ray rates for the whole field of view. The fourth line (drawn at 10^ 
times its true value) gives the number of cosmic ray events (after gamma cuts 
and including energy bias) within the source region, which is defined by the 
angular resolution shown in figureU(using no jitter). 



can be used 114411 . Usually, the background is measured in a sky 
region within the same field of view, but sufficiently far away 
from the source to avoid a contribution from the source itself. 
This background region can be larger than the source region to 
reduce the statistical uncertainty on the background estimate. 
The ratio of the solid angles of the source and the background 
regions is expressed by the a-factor: 



(10) 



Note, that the choice of the optimum size of the source region depends on 
the source strength as well as its energy spectrum which are a priori not known. 



^background 

For the conservative estimate, a = 1 is assumed. For the 
optimistic scenario, a <K 1 is assumed, which means that the 
background region is chosen much larger than the source re- 
gion. Due to the large field of view of HiSCORE, and the rather 
flat instrument acceptance (see 12811 for acceptance plots), a low 
a-factor is anticipated. 

Using the background levels shown in figure [12] the number 
of gamma-rays with energy E > Eq needed to fulfil the detec- 
tion criteria is calculated for a range of values of Eq between 
50 TeV and lOPeV. The gamma-ray flux is modelled using a 
spectrum proportional to E^^-^ without a cutoff, and the flux 
constant is adjusted to produce the required number of gamma- 
rays. 

Figure [13] shows the calculated sensitivity for the HiSCORE 
detector in comparison with other planned or existing gamma- 
ray observatories. For HiSCORE, three scenarios are plotted: 
A lOkm^ array with conservative assumptions (1 ns time res- 
olution, a - 1), a lOkm^ array with optimistic assumptions 
(ideal time synchronisation, a <k 1), and a lOOkm^ array with 
optimistic assumptions. Due to the wide field of view of 0.6 sr, 
about 25% of the sky will receive an exposure of more than 200 
hours per year [23]. Therefore, the sensitivity has been calcu- 
lated for an observation time of 1000 hours, equivalent to five 
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Figure 13: Point source sensitivities for tlie simulated HiSCORE detector 
("HS"), for conservative and optimistic assumptions (see section l5?2l for de- 
tails). Also shown is the expected sensitivity for a 100 km- HiSCORE array. 
For comparison, published sensitivities of selected other gamma-ray observa- 
tories are shown (CTA |45], H.E.S.S. |46], Milagro and HAWC |47]). The 
observation time is assumed to be 50 hours for pointed instmments (H.E.S.S. 
and CTA) and five years of continuos operation for all other instruments. 



years of continuous operation. 

It can be seen that HiSCORE can contribute sensitive obser- 
vations in a so far poorly covered energy (or wavelength) band 
of the electromagnetic spectrum. A more detailed discussion 
of the scientific potential of an instrument with this sensitivity, 
also in the light of existing measurements and models, can be 
found in iH]. 



6. Conclusions 

An event reconstruction algorithm for a non-imaging air 
Cherenkov array with sparse sampling has been presented, and 
its performance has been evaluated using a dedicated simulation 
framework. The choice of station spacing and detector charac- 
teristics have been guided by the HiSCORE array, but the meth- 
ods should be applicable as well to other air shower detectors 
of similar kind. The algorithm incorporates techniques well es- 
tablished in the field, but also new approaches like the arrival 
time model for the direction and shower depth reconstruction 
or the signal stacking for the measurement of the signal width. 
The gamma hadron separation algorithm has been newly de- 
veloped, partly using ideas that have been raised in other stud- 
ies. The described algorithm has been used to demonstrate that 
even with a sparse array such as the simulated one (station spac- 
ing 150 m) an accurate event reconstruction is possible down to 
about 50 TeV (see table©. 

The results show that with the assumed configuration the 
HiSCORE detector can be a powerful observatory for gamma- 
ray astronomy. Its sensitivity is sufficient to study the continua- 
tion of currently known gamma-ray source spectra to ultra high 
energies, and offers a great discovery potential for Galactic pe- 
vatrons and other sources of ultra high energy gamma radiation. 





at 


100 TeV 


at 1 PeV 


Core position res. [m] 




20 


5 


Direction res. [deg] 




0.25 


0.1 


Energy res. [percent] 




20 


10 


Depth res. [g/cm^] 




70 


40 


Sensitivity lOkm^ [ergs^'cm^^] 


4 


X 10-'^ 


5 X 10-'^ 


Sensitivity lOOkm^ [ergs^'cm^^] 


1 


X 10-'^ 


6x lO-i'* 



Table 1: Performance figures of the described reconstruction algorithm, as 
derived with the HiSCORE simulations. Core, direction and energy resolu- 
tions are given for gamma-ray events, depth resolution for proton events. The 
HiSCORE sensitivities are given for the optimistic scenario. 
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